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Abstract 

We formulate a self-consistent field theory for polyelectrolyte brushes in the presence 
of counterions. We numerically solve the self-consistent field equations and study the 
monomer density profile, the distribution of counterions, and the total charge distri- 
bution. We study the scaling relations for the brush height and compare them to the 
prediction of other theories. We find a weak dependence of the brush height on the graft- 
ing density. We fit the counterion distribution outside the brush by the Gouy-Chapman 
solution for a virtual charged wall. We calculate the amount of counterions outside the 
brush and find that it saturates as the charge of the polyelectrolytes increases. 

1 Introduction 

Colloidal particles have various industrial applications. Their tendency to flocculate some- 
times causes problem in applications. Charging colloidal particles is an effective way to 
prevent the floccuiation by electrostatic repulsions. However, such charge effects are easily 
screened by salt and polar solvents in the surrounding environment. Grafting polymers onto 
the colloid surface offers another way to prevent the floccuiation. Steric repulsion induced by 
the grafted polymers makes colloidal particles repel each other. 

Grafting polyelectrolytes onto the colloid surface offers an even more efficient way to pre- 
vent the floccuiation of the colloids in polar solutions [I]. The stabilization effect arises from 
both steric and electrostatic repulsion. Since the polyelectrolytes can trap their own coun- 
terions in their vicinity, the grafting of polyelectrolytes generates a layer of locally enhanced 
counterion concentration. The stabilization effect obtained by grafting polyelectrolytes is thus 
less sensitive to the salinity of the surrounding environment than the one obtained by solely 
charging the colloids. 
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Grafted neutral polymers have been studied by Alexander [2] and de Gennes[5]. When 
polymer chains densely grafted onto an impenetrable planar surface are immersed in a good 
solvent, the excluded volume effects strongly stretch the chains and they form a blob structure. 
This particular distribution of densely grafted polymers is called a "brush". In this neutral 
brush regime, they found that the density profile of the brush is essentially flat and they 
predicted a linear dependence of the brush height on the polymerization index and some 
other scaling for the grafting density. 

Later, self-consistent field theories were developed for neutral polymer brushes [1]- [8]. In 
these papers, it was shown that polymer chains are not uniformly stretched, and the density 
profile of the brush is not flat but parabolic. 

Polyelectrolyte brushes which consist of charged polymer chains densely grafted onto a sur- 
face have been extensively studied by theories [H]-[23], simulations [24|-[3"U] and experiment s [3 1 j - 
|48j . The free energy of the polyelectrolyte brush consists of three terms in general: The 
osmotic pressure of counterions, the Coulomb interactions between the charges, and the ex- 
cluded volume interactions of monomers. 

Pincus [IT] has proposed a scaling theory for the polyelectrolyte brush assuming that its 
equilibrium height is determined by balancing the osmotic pressure of the counterions with 
the Gaussian elasticity of the chain. Scaling laws were presented for two regimes: (1) the 
osmotic regime: When the polyelectrolyte brushes are densely grafted and have a large value 
of their charge fraction (strongly charged case), most of the counterions are captured inside 
and in the vicinity of the brush due to the strong Coulomb attraction. The brush height is 
roughly equal to the thickness of the counterion layer and it is independent of the grafting 
density. (2) Pincus regime: When the chains are sparsely grafted and have small values of 
their charge fraction (weakly charged case), the Coulomb attraction is relatively weak and 
the counterions are distributed far beyond the brush region. The counterion swelling pressure 
is reduced by the fraction captured inside the brush. 

The polyelectrolyte brush has been further investigated with the box model where a flat 
density distribution is assumed, and a collapsed regime was introduced [2TJ [28]. In this 
regime, attractive forces between grafted polyelectrolytes occur due to the correlations of the 
monomers and counterions. The brush height is determined by balancing the above attractive 
forces and the repulsive forces due to the excluded volume effects. 

Zhulina et al. [2] developed a self-consistent analytical model and considered the effect 
of the excluded volume interactions in the osmotic regime. In the osmotic regime, for good 
solvents, they obtained a scaling relation which does not agree with that of the box model. 
This work was generalized to poor solvents [131 CHI- Borisov et al. [T7j introduced a quasi- 
neutral or Alexander regime where the excluded volume effects dominate over electrostatic 
effects. The quasi-neutral regime has the same scaling relation as the neutral brush except 
for the dependence on the charge fraction of the chain. The brush height is determined by a 
balance between the excluded volume effects and the chain elasticity. 

Self-consistent field theories have been developed for polyelectrolytes by several authors. 
Borukhov et al. applied a self-consistent field theory to semidilute solutions of polyelectrolytes 
and polyampholytes [H]. Wang et al. applied it to the interface of phase-separated polyelec- 
trolyte solutions E3]. Witte et al. applied it to the mixed polyelectrolyte brushes [2"3]. 

The phase diagram for various values of the charge fractions and grafting densities has 
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been obtained from scaling theories [T71 |2"T] or simulations |29]. The Bjerrum length Zb is a 
characteristic length scale at which the electrostatic interactions between monovalent ions and 
the thermal energy become equal. When Z B > v, where v is the excluded volume parameter, 
the electrostatic interactions are relatively strong. This situation is called the 'strong coupling 
case'. When Z| < v, the electrostatic interaction is relatively weak. This situation is called 
the 'weak coupling case'. Scaling theories [21] and molecular dynamics simulations [28] both 
predict that a collapsed regime appears in the strong coupling case. 

In the present work we formulate the self-consistent field theory for a salt-free poly- 
electrolyte brush in the presence of counterions based on a theory developed by Orland 
and Schick [6] for the neutral polymer case combined with a theory developed by Netz and 
Orland [SU] for charged electrolytes. We apply it to polyelectrolyte brushes grafted onto an 
impenetrable planar surface. The chain configurations are obtained by solving some mod- 
ified Edwards equations, including electrostatic interactions. The electrostatic potential is 
obtained by solving the nonlinear Poisson-Boltzmann equation. 

We solve these self-consistent field equations numerically in the weak coupling case and 
obtain the monomer density profile, counterion distribution, and total charge distribution. 
These results are obtained for various charge fractions and we study the scaling relations 
resulting from our model. 



2 The Partition Function for Polyelectrolyte Brushes 



We consider M charged polymers of polymerization index N grafted to a planar surface of 
area A. We assume that the grafting surface is impenetrable so that all polymers are confined 
to a half-space. The monomer positions are denoted by {Rfc(s)}, where s G [0, N] is the 
curvilinear abscissa along a chain of N monomers and k = 1, • • • , M is the label of the M 
chains. In addition, we assume that the s = extremities of the chains are grafted to the 
surface, at positions Rfc(0) = (u k ,z k (0) = 0), where z denotes the axis perpendicular to the 
grafting plane. We consider N c counterions and denote their positions rj, where i — 1, • • • , N c . 
The canonical partition function Z^ c for this system can be written as 



J N C 
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where / T>g is the functional integral over the function g(r). The first term in the exponent 
represents the elasticity of the Gaussian chain, where Rjt(s) is the derivative of Rfc(s) with 
respect to the curvilinear abscissa s and a is the Kuhn length. The grafting of the chains is 
imposed by the boundary conditions. 

The second term represents the electrostatic interaction where k B T is the thermal energy, 
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is the Coulomb interaction, e is the dielectric constant of the solvent and p c (r) is the local 
charge density including all charges in this system (charged monomers and counterions) : 



M „ N N c 

^2 ds q k (s)5(r - R k (s)) + ^ ge5(r - r 
k=i Jo i=i 



(3) 



where qk{s) is the variable denoting the charge carried by the monomer s along the chain k. 
The charge of the counterions is denoted by qe. Without loss of generality we assume that all 
polymers are negatively charged. We use the smeared model for the charge of the chains. For 
a polyelectrolyte with a fraction / of its monomers negatively charged, this model assumes 
that each monomer carries a uniform fractional charge — fe, where — e is the electron charge. 
Namely, qk{s) = —fe for any monomer s on any chain k. The condition of electroneutrality 
is then expressed as 

J dr Pc (r) = 0, (4) 

which yields N c = fMN. 

The third term in equation ([1]) represents the excluded volume interactions between 
monomers, where v is the excluded volume parameter and 



M rN 



An(r) = J2 ds5(r- Rfe(s)) 
k=i Jo 



(5) 



is the local monomer concentration. 

The operator inverse of the Coulomb interaction for a system with homogeneous dielectric 
constant e is given by 

^(r,!- 7 ) = -eV 2 5(r-r'). (6) 
We introduce two fields <pi, <p2 to perform the Gaussian transformations 
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and the partition function reads 
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where l B = e 2 /4irek B T is the Bjerrum length. 
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The partition function is brought into a more tractable form by going to the grand canon- 
ical ensemble 



(10) 



N c =0 



where A is the fugacity which is related to the particle chemical potential fi by A = e M//fcsT . 
The grand canonical partition function can be written as 



VfaVfaexp [-H [0!,0 2 ]] 



with the Hamiltonian 

H[<f>xM = 

where 

Qk [4>i, 4>2] 
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k=l 



R(0)=(u fc ,0) 



DR(s) exp 



A 1 



ds { t^R 2 (s) - t-MK(s)) + i^(R(s)) 



2a 



q 



(13) 

As follows from the definition of the grand canonical partition function, equation ffTOl) . the 
expectation value of the number of counterions is given by 



(N c ) = A 



d\nZ> 



A( / dre- iMr) ). 



(14) 



In the limit M — > oo, if the grafting points on the surface are random independent 
variables, we can further simplify eq. fll2l) . Indeed, according to the central limit theorem, we 
have 



^lnQ fc [0i, (j) 2 \ ~jvf— > 
fc=i 



M J rfuP(u) lnQ[0 1 ,0 2 ] + O(v / M), 
where -P(u) is the probability distribution of the grafting point u and 
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In the following, we will assume a uniform grafting density, that is 

1 



(17) 



where A is the total grafting area. If we define a = A/M as the average grafting area per 
chain, we have 



V(piV(f)2 exp 



dr 



[V0i(r)] : 
%nl B q 2 



+ -^ 2 2 (r) - Ae~^ r A + 1 JjulnQifa, 



(U 



5 



3 Mean-field Approximation 



In order to evaluate the partition function (|T8|) . we use the saddle-point method to determine 
the fields <p\ and 02 which most contribute to the functional integral. We assume that these 
two fields possess translational symmetry in the (x, y) directions, parallel to the grafted plane 
(remember that z is the direction normal to the plane). The Hamiltonian reads 



H[(t>u 



A / dz 



Ml B q 7 



' dcj)i(z) 
dz 



(19) 



where 



Q 



z(0)=0 



Vz(s) exp 



N f 3 / 

ds | 7^z 2 ( s ) - *-0iOO)) + 20 2 (z(s)) 



(20) 



The saddle-point method provides the mean-field approximation to the system. The mean- 
field $2(2) are given by the stationarity conditions 



1> V2 



<tyl(*) 
^[01, 2 ] 



0. 



0. 



(21) 
(22) 

(23) 



The mean-field approximation to the Hamiltonian H m f is simply 

H mi = H[i<$> 1: i<S> 2 ]. 

Note that $1(2) is identified as the electrostatic potential. 
The mean-field equations are thus 

1 d^jz) + Xe ^ i{z) _ 1 5\nQ = Q 
47r/ B g 2 cfe 2 a5$i(z) 

v 5\nQ v 

where a = A/M is the average area per chain and p m (z) is the local density of monomers. 
The density of monomers can be written as 

lr r N r r N ( 3 f t> 1 

Pm(^) = - Vz(s) ds5(z-z(s))exp - ds I —z 2 (s) - ^^(s)) + -p m (z(s)) \ . 
Q J Jo Jo { 2er q a J 

(26) 

The density is normalized so that 



(24) 



(25) 



p m {z)dz = N. 



(27) 
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We take the grafting plane to be at z = and introduce the monomer probability ampli 
tudes 

1>(z,s) = {z\e-' H \0), (28 
tp^(z,s) = 



dz'(z'\e- sH \z) 

where we use the standard quantum-mechanical notations 
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• (30) 



Note that as before in the paper, the curvilinear abscissa s is defined so that its values at the 
grafted and free ends are and N respectively. 

The probability amplitudes if)(z, s), ip^(z, s) satisfy the modified diffusion equations 
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ip(z,s) = 0, 
?p\z,s) = 0. 



(31) 
(32) 



with the initial conditions ip(z,0) = S(s) and ip^(z,0) = 1. Because the grafting plane is 
impenetrable,-*/^, s), ip^(z, s) must satisfy the boundary conditions -^(0, s) = 0, ^(0, s) = 0. 
The density of monomers is given by 



Pm(z) 



J N dsi/j(z,s)^(z,N 



and equation (124)) can be written as 

1 d 2 $i(z) 



Anl B q 2 dz 2 
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(33) 



(34) 



This is just the Poisson-Boltzmann equation for all the charges in the system (polymers and 
counterions) . In the following, we will consider only monovalent counterions, that is q = 1. 
The counterion density per unit area of the grafting plane is given by 



p ci (z) = \e-*^\ 



(35) 



From equation (134")) one obtains / °° p c \(z)dz = (N)/A. The condition of electroneutrality is 
given by 

fN 

p c i(z)dz = . (36) 



a 



Equations ( 1311 to (1361) are the self-consistent field equations describing the polyelectrolyte 
brush in the presence of counterions. We assume that the grafting plane is electrically neutral. 
The electroneutrality of the system is enforced by setting the field strength to zero at the 
boundaries 



dz 



z=0 



d^t(z) 
dz 



0. 



(37) 



The total charge density p c (z) per unit area of the grafting plane is given by 

f 

Pc(z) = PmO) + PciO)- (38) 

a 

The total charge distribution p c (z) satisfies the condition of the electroneutrality 

POO 

/ p c (z)dz = 0. (39) 
Jo 

4 Numerical Solutions of the Self-Consistent Field Equa- 
tions and Discussions 

In this section we present numerical solutions of the self-consistent field equations as a function 
of the values of the polymerization index N, grafting area per chain a, excluded volume 
parameter v, charge fraction /, and Bjerrum length In this paper, we vary / from to 
1.0, while we fix the values of parameters iV = 30, v = 1.0a 3 , Zb = 0.1a, a = 10a 2 , otherwise 
stated. This corresponds to the weak coupling case Z B < v where the Coulomb interaction is 
relatively weak. To investigate the grafting density dependence, we vary a from 5a 2 to 20a 2 . 
To investigate the polymerization index dependence, we vary iV from 30 to 70. 



4.1 Structure of the Brush 

In Figured] we show the monomer density profiles p m {z) for the above fixed values and various 
values of / including the monomer density profile of the neutral brush (/ = 0). The brush 
height increases with increasing /. 

The impenetrability of the grafting plane induces a depletion area near the grafting surface 
in the monomer density profiles and pushes the maximum of the density away from the 
grafting surface. The depletion area is also observed for the neutral brush[5T]. The density 
maximum gets closer to the grafting plane with increasing /. This is due to the strong 
Coulomb interaction between the monomers and the counterions for large /. The monomer 
density decreases exponentially when the distance from the grafting plane z becomes large. 

Taniguchi et al. [52] obtained similar results but the behavior of the monomer density near 
the grafting plane is somewhat different from our results. 

We measure the brush height h m and the thickness of the counterion layers l„ as the first 
moments of each density profiles: 

fjfzpj^dz 
/ °° zp ci (z)dz 

cl roc f \ j • v^ 1 / 

j p ci {z)dz 

We show the values of 2h m and 2Z C i for each values of / in Table [U We multiply by a factor 
2 because the height for the uniform density distribution corresponds to twice of the first 
moment. 
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In Figure [2k we plot h m and / on a log-log scale. For large values of / > 0.5, the brush 
height behaves as h m ~ f a , with exponent a = 0.41. In Figure [2b we plot h m and 1 + / on a 
log-log scale. Figure [2b shows the scaling relation h m ~ (1 + f) 13 with exponent /? = 1.17 for 
small values of / < 0.5. 

We investigate the dependence of the brush height h m on the other parameters N and a. 
Figure [3] shows the scaling relation h m ~ iV 7 with exponent 7 = 1.0 for / = 0.3, 0.5, 1.0 respec- 
tively. Figure [4] shows the scaling relation h m ~ (1/cr) 5 with exponent 5 = 0.06, 0.13, 0.18, 0.23 
for / = 1.0,0.5,0.3,0.0 respectively. 

We briefly review the scaling relation predicted by the scaling theory and other SCF 
theory. For the strongly charged case, Pincus[TT] predicted that the brush height h behaves 
as 

h ~ Naf 1/2 , (42) 
a = 1/2,7 = M = 0, 

in the osmotic regime. 

Zhulina et al. [H] found that the scaling relation in the osmotic regime changes due to 
excluded volume effects to 

h ~ Naf 2/5 , (43) 
a = 2/5,7 = 1,5 = 0. 

The scaling relations ( 1421) . (1431) indicate that the brush height h is independent of the grafting 
density 1/er. 

Recent study for the osmotic regime shows that the brush height h slightly increases with 
increasing l/q[22l [29l 1471 148] . The term 'nonlinear osmotic regime' is used for this weak 
dependence on 1/er. (5 = 1/5 for weak coupling and / = 1) 

PincuspT] also obtained the brush height in the Pincus regime where counterions are 
distributed beyond the brush behaves as 

h ~ N 3 a 2 f 2 /a, (44) 
a — 2, 7 = 3, S — 1. 

The quasi-neutral regime [21 [31 [17], [21] appears both for the strongly and weakly charged 
case. The brush height is determined by balancing the repulsive force due to the excluded 
volume effect and the attractive force due to the chain elasticity. The scaling relation is 
identical to the neutral brush case except for the dependence on / for the strongly charged 
case. The scaling relation for the strongly charged case is 

h ~ (1 + f) 2/3 N(a 2 v/a) 1/3 , (45) 
(3 = 2/3,7 = M =1/3. 

For the weakly charged case, the scaling relation is exactly the same as that for the neutral 
brush: 

h ~ N^v/a) 1 / 3 , (46) 
a = 0,7 = 1,5 = 1/3. 
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We compare our numerical results with the other theoretical predictions. First we inves- 
tigate the results for / > 0.5. Our numerical results a = 0.41 is in good agreement with 
a = 2/5 of (j43|) for the osmotic regime in the good solvent. Our numerical results 7 = 1.0 
for / = 0.5, 1.0 is consistent with (I43p . The brush height scales linearly with N. Our results 
show that the brush height h m depends weakly on the grafting density 1/er with an exponent 
5 = 0.06,0.13 for / = 1.0,0.5 respectively. The scaling relation (14"3"|) predicts 5 = 0. These 
comparisons show that our numerical results for / > 0.5 corresponds to the osmotic regime. 
The slight increase of the brush height for increasing the grafting density is consistent with a 
recent study of the nonlinear osmotic regime. [22j [29J, HH HH] . 

Next we investigate our results for / < 0.5. Our results of (3 = 1.17 does not agree with 
the value f3 = 2/3 of the scaling relation (1461) for the quasi- neutral regime. Our results of 
7 = 1.0 for / = 0.3 is consistent with (JMP . Our results of 5 = 0.18, 0.23 for / = 0.3, 0.0 differ 
from the exponent 5 = 1/3 of (]46j) . 

The neutral brush is characterized by a single stretching parameter [3, E] 



For large value of b > 10, the monomer density profile is well approximated by the semi 
classical parabola density profile [5] and the brush height satisfies the scaling relation (147j) . 
However the deviations from the parabolic profile become important when the values of 
N,l/(x,v decrease and the stretching parameter becomes b < 10. Our values of the pa- 
rameters N = 30, v = 1.0a 3 , 5a 2 < o < 20a 2 correspond to 11.7 > b > 4.7. Therefore, the 
disagreement between the calculated value of the exponent 5 (0.06, 0.13, 0.18, 0.23) and the 
value of the scaling prediction ( T4T1) (8 = 1/3) is not surprising. 

When / decreases, the value of the exponent S becomes close to the value for the neutral 
case. Roughly speaking our results for small value of / < 0.3 correspond to the quasi-neutral 
regime of the phase diagram [TT], ETJ [29] . 

Our value for the exponent (3 = 1.17 does not agree with the value (3 = 2/3 of the scaling 
relation (1461) . When (146j) was derived, only excluded volume interactions were considered 
and the osmotic pressure of counterions was neglected. A more realistic value of the expo- 
nent (3 should be derived when the excluded volume interaction and the osmotic pressure of 
counterions are equally treated. 

Our numerical results show a crossover phenomena between the osmotic regime and the 
quasi-neutral regime. Scaling theories assumes that one of the two repulsive effects is dominant 
and balance it to the elasticity of the Gaussian chain. But in many regimes the two effects 
are comparable; The behavior of the brush does not change abruptly and a smooth transition 
between the two regime is expected. Our numerical results for 0.4 < / < 0.5 correspond to 
this crossover region. 

4.2 Distribution of Counterions and Charges 

Next we consider the distribution of counterions and charges. We first plot the distribution 
of counterions p c \(z) for various values of / in Figure [5j We observe a depletion area of Pa{z) 
near the grafting surface. This depletion effect is clear for / > 0.3. 




(47) 
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At large distance from the wall, the monomer density decreases faster than the counterion 
density. Therefore there are no more monomers in the solution, and the counterions see a 
virtual charged wall with an effective charge density a e s at position r D g. 

We study the thickness of the counterion layer. In Table [fl we show the values of 2Z C i for 
each values of /. 

We fit the counterion distribution at large distance by the GC profile 

/(*) = tttt ni . (48) 
2tiIb{z — Li) 1 

where Lf is a fit parameter. 

For z ^> Lf, eq. (!48p is equivalent to the distribution of counterions in a system which 
consists of a virtual charged plane located at r cff with charge density — a^e in the presence 
of monovalent counterions. The parameter r e g is set to be the minimum value which satisfies 
the condition „ r ff 

" p c (z)dz = - / f(z)dz, (49) 

which expresses the fact that the total charge to the left of r e g is equal to that of the coun- 
terions on the right of r e ff. cr e ff is determined by 

acS = o / 1 1 ETT' ( 50 ) 

27r/ B (r cff - H m ) 

The values of r e s, cr c s are presented in Table [TJ The effective charge density a c ^e is almost 
equal to the integrated amount of charges of counterions for z > r e g. The total amount of 
charges of counterions is a s = fN/a. From a e g we define the effective Gouy-Chapman (GC) 
length as Z e ff = 1/27tIb<J c s- This effective GC length Z e g provides a definition of the thickness 
of the counterion layer outside the virtual wall. We may thus take L ci = r e g + l c g as an 
alternative definition of the thickness of the counterion layer. 

In Figure [6] and Figure [TJ we plot the absolute value of the negative charge distribution 
fN/ap m (z), the positive charge distribution p c i(z), the total charge distribution p c (z), and 
the fitting GC solution (T4"g|) . In Figure [H] and O we plot cr c{i /a s , L ci /r cS . 

For large values of / > 0.5, the ratio cr cS /a s is smaller than 0.10 and a large part of the 
negative charge of the brush is cancelled by the positive charge of counterions. This strong 
screening of the negative charges of the brushes can be verified by comparing Figure [10] and 
[6l For decreasing /, the screening of the brush charge becomes weaker and the counterions 
gradually escape from the brush. This behavior of counterions is verified in Figure [71 

From Tabled], we see that for / > 0.2, the value of <7 e g ~ 0.12 is almost constant. As was 
already mentioned above, far away from the wall, the polyelectrolyte brush can be viewed 
as a charged plane located at r e g with charge density — a e ge. The charge density — a c ge is 
the 'effective' charge of the brush which is not cancelled by the counterions. Our numerical 
results show that there is an upper limit for this effective charge. 

Next, in Figure [10] we plot the total charge distributions p c (z) for various values of /. 
Figure [TU] shows that the values of p c {z) near the grafting plane are always positive for any 
/. This is due to the depletion of the monomer density near the grafting plane. The amount 
of positive charges of counterions is larger than that of the negative charges of the brushes. 
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The minimum of p c (z), which is negative, is located near the grafting plane and an electric 
dipole layer is generated. 

The depth of the negative minimum of p c (z) increases with / and a small peak in the 
absolute value of the positive charge (fN/a)p m (z) appears for / = 1.0 (see Figure [6]). As 
can be seen in Figure El when / increases, in order to enforce locally charge neutrality, the 
monomer and counterion densities become closer. However, the depletion layer of the polymer 
near the wall is still present. This region is filled by the counterions, and this accumulation 
of positive charges is compensated by a peak in the polymer density. Taniguchi et al. [52] 
obtained similar results. However the behavior of the depth of their negative minimum is 
somewhat different from our results. 

5 Conclusion 

We formulated the self-consistent field theory for a planar polyelectrolyte brush in a salt free 
solution in the presence of counterions. We numerically solved the self-consistent field equa- 
tions for the weak coupling case and obtained the monomer density profile, the distribution 
of counterions, and the total charge distribution. We also studied the scaling relation for the 
brush height. For large values of the charge fraction, our results correspond to those of the 
osmotic regime. We observe a weak dependence of the brush height on the grafting density. 
This is consistent with a recent study of the nonlinear osmotic regime. For small values of 
the charge fraction, our results correspond to those of the quasi-neutral regime. We find a 
crossover phase between the osmotic regime and the quasi-neutral regime. 

We studied the distribution of the total charge, and found that an electric dipole layer 
appears near the grafting area. We studied the distribution of counterions, and fit it by the 
Gouy-Chapman theory for a virtual charged wall at the location r e g with an effective charge 
density a c g . We calculated the integrated amount of the charge of counterions outside the 
brush and the thickness of the counterions outside the brush by introducing an effective GC 
length (leg). For large values of the charge fraction, the charges of the monomers are strongly 
screened by the charge of the counterions, and the counterions are trapped inside the brush. 
For smaller values of the charge fraction, the screening effects are weak and the counterions 
extend far beyond the brush. Far away from the grafted surface, the polyelectrolyte brush can 
be viewed as a virtual charged plane with an effective charge density a e g. We found that this 
effective charge density has an upper limit and the rest of the polymer charges are cancelled 
by the charges of the counterions. 
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Table 1: 2h m , 2l c i, r e g, l e g, L c ;, <t s , a c g at different charge fractions /. The other parameters are 



the same as in Figured! 



/ 


2h m /a 


2/d/a 


W a 


ks/a 


L C i/ a 


0"s 


Ocff 


1.0 


18.61 


26.16 


30.96 


14.22 


45.18 


3.0 


0.112 


0.9 


17.78 


26.12 


29.26 


13.98 


43.24 


2.7 


0.114 


0.8 


16.91 


26.26 


27.49 


13.74 


41.23 


2.4 


0.116 


0.7 


16.00 


26.62 


25.67 


13.49 


39.16 


2.1 


0.118 


0.6 


15.05 


27.36 


23.84 


13.31 


37.15 


1.8 


0.120 


0.5 


14.05 


28.70 


21.96 


13.17 


35.12 


1.5 


0.121 


0.4 


13.01 


31.10 


20.01 


13.09 


33.10 


1.2 


0.122 


0.3 


11.92 


35.64 


18.00 


13.19 


31.19 


0.9 


0.121 


0.2 


10.83 


45.36 


16.18 


13.97 


30.15 


0.6 


0.114 


0.1 


9.81 


74.52 


14.41 


16.76 


31.17 


0.3 


0.0950 
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Figure 1: Monomer density profile Pm(z) for polymerization index N = 30, grafting area per 
chain a = 10a 2 , excluded volume parameter v = 1.0a 3 , Bjerrum length 1-q = 0.1a. Profiles for 
the charged fraction /=0 (neutral), 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 1.0 are shown. 
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Figure 2: (a) The brush height h m as a function of the charge fraction /. The parameters are 
the same as in Figured! We plot / = 0.1, 0.2, •••,(). 9, 1.0 from left to right, solid line shows a 
power law fit. (b)The brush height h m as a function of 1 + /. The parameters are the same as 
in Figure [TJ We plot / = 0.1,0.2,---,0.9,1.0 from left to right. The solid line shows a power 
law fit. 
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Figure 3: Brush height h m as a function of the polymerization index N. f = l.O(circles), 
/ = 0.5 (triangles), / = 0.3(squares) are shown. The other parameters are the same as in 
Figure [TJ 
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Figure 4: Brush height h m as a function of the grafting density I /a. f = l.O(circles), 
/ = 0.5 (triangles), / = 0.3(squares), / = O.O(cross). The other parameters are the same as 
in Figure [TJ 
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Figure 5: The counterion distribution p&(z) for / = 0.1,0.2,0.3,0.4,0.5,0.6,0.8,1.0. The 
other parameters are the same as in Figure [H 
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Figure 6: The thick solid lines are absolute values of the negative charge distributions 
fN/ap m (z), the thick dashed lines represent the positive charge distribution Pd(z), the thick 
dotted lines are the total charge distributions p c (z) , and solid lines are fitting curves for the 
Gouy-Chapman solution. The parameters are the same in Figure [TJ (a) f=1.0, (b) f=0.6, 
(c)f=0.4. 
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Figure 7: The thick solid lines are absolute values of the negative charge distributions 
fN/ap m (z), thick dashed lines are the positive charge distributions Pd(z), the thick dot- 
ted lines are the total charge distributions p c (z) , and the solid lines are fitting curves for 
the Gouy-Chapman solution. The parameters are the same in Figured! (a) f=0.3, (b) f=0.2, 
(c)f=0.1. 
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Figure 8: cr c ff/er s for / = 0.1, 0.2, • • • , 0.9, 1.0. The other parameters are the same as in Figure 
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Figure 9: L C i/r c g for / = 0.1, 0.2, • • • , 0.9, 1.0. The other parameters are the same as in Figure 

m 
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